Binding of holons and spinons in the one-dimensional anisotropic t-J model 
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We study the binding of a holon and a spinon in the one-dimensional anisotropic t-J model using a Bethe- 
Salpeter equation approach, exact diagonalization, and density matrix renormalization group methods on chains 
of up to 128 sites. We find that holon-spinon binding changes dramatically as a function of anisotropy parameter 
a = J±/J z : it evolves from an exactly deducible impurity-like result in the Ising limit to an exponentially 
shallow bound state near the isotropic case. A remarkable agreement between the theory and numerical results 
suggests that such a change is controlled by the corresponding evolution of the spinon energy spectrum. 
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In the one-dimensional t-J and Hubbard models, spin and 
charge dynamics are independent, leading to the well-known 
effect of spin-charge separation: the splitting of the elec- 
tron (hole) into spinon and holon excitations that carry only 
spin and only charge, respectively flj]. Considerably less is 
known about interactions among these excitations. Recently, 
it was shown that in the supersymmetric t-J model, spinons 
and holons attract each other [2], but this does not result in 
a bound state, indicating that the spinon-holon interaction is 
non-trivial. 

On the other hand, there is a limit in which spinon and holon 
do form abound state: the t-J z model jj. Here, the half-filled 
system is an Ising antiferromagnet (AF). Removing a spin and 
moving the hole away from origin increases the net magnetic 
energy by J z /2 due to creation of a domain wall in the AF or- 
der - an immobile spinon (Fig. [T). Once moved, the hole also 
carries an AF domain wall and is a free holon with dispersion 
el = — 2icosfc. Since recombination with the spinon lowers 
the energy of the system, the holon can be considered as mov- 
ing freely except at the origin where it is subject to an effective 
attractive potential F = — J z /2. Clearly, such a potential in ID 
always leads to a bound state. The single-hole Green's func- 
tion can be found exactly: G{w) = [J z /2 ± y/u> 2 — 4t 2 \ 
Its lowest pole gives the binding energy of the spinon-holon 
bound state: 
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in which holon is confined in the vicinity of the spinon with 
the localization length t oc A -1 / 2 . This leads to the question: 
how does the Ising limit, in which the spin and charge are 
bound, evolve into the isotropic limit where they separate? 
This question, together with the expectation that the study of 
anisotropic model can shed light on the general problem of 
spinon-holon interactions, motivates this work. 
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FIG. 1: (Color online). A hole in the Ising AF (circle) moved by four 
sites from origin. The spinon is indicated by the dashed box. 



In this Letter we address this problem by a detailed theoreti- 
cal and numerical study of the anisotropic t-J model using the 
Bethe-Salpeter equation (BSE), exact diagonalization (ED), 
and density matrix renormalization group (DMRG) [4] meth- 
ods on systems of up to 128 sites. Although the binding be- 
comes too weak to be detected numerically near the isotropic 
limit, we believe that a holon and a spinon form a bound state 
for any value of anisotropy a = Jj_/ J z less than one. Given 
a remarkable agreement between the theory and numerical re- 
sults, we argue that the binding is largely controlled by the 
spinon energy spectrum. Qualitatively, at a <C 1 the holon 
is confined in the vicinity of the impurity-like spinon, while 
already for a > 0.5 the spinon is quasi-relativistic and is only 
weakly bound to the holon. Such a change in the pairing char- 
acter results in a dramatic decrease of the binding energy at 
larger a. In addition, using a small-a analysis, we provide 
insight into the kinematic structure of the spinon-holon inter- 
action that offers a simple explanation for having an attraction 
but no bound state in the isotropic limit of the model. 

Using standard notation, we define a generalization of the 
ID t-J model, H 
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where i and j = i ± 1 are the neighboring sites. This model 
interpolates between the t-J z model (a = 0) and the isotropic 
t-J model (a = 1). We set t = 1. 

To obtain the binding energy numerically, we have per- 
formed calculations of the ground state energies (GSEs) us- 
ing both ED and DMRG techniques. While only small system 
sizes are accessible by ED (up to L = 23 sites in our case), 
it allows us to calculate the lowest energy in each momen- 
tum sector (Fig. [2]) and is also important in analyzing various 
finite-size effects. Using DMRG [4] we have calculated GSEs 
of systems of up to L — 128 sites using periodic boundary 
conditions (PBCs), which greatly increases the numerical ef- 
fort required. Up to m = 1400 states per block were kept 
in the finite system method, with corrections applied to the 
density matrix to accelerate convergence with PBCs |5|]. The 
binding energy A for the infinite system was obtained using 
two different extrapolation methods described below. In the 
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FIG. 2: Lowest energies Ek—Eo vs k in an L = 21 chain with (a) zero 
and (b) one hole, J z =4.0. Solid lines in (a) is the spinon dispersion 
from BA [9]; in (b) the lines are guides to the eye. 




FIG. 3: (Color online). Theoretical (lines) and numerical results by 
method A (circles) and method B (diamonds) for A vs a, J z = 4.0. 
Inset: | A| for a > 0.4 on a semi-log plot. 



range of a where the calculation of A is reliable (a ^0.6), we 
found excellent agreement between the numerical data from 
both methods, and theoretical results, based on BSE. A sum- 
mary of the results for J z = 4.0 is presented in Fig. |3j results 
for other values of J z are qualitatively similar. 

Numerics. For a given J z and a, the binding energy is 



A — E p — E s — Eh, 



(3) 



where E x , x — p,s,h, is the GSE of the L = oo chain with 
the spinon-holon pair, the spinon, and the holon, all relative 
to GSE of the XXZ chain, respectively. It is well known 
10,13] that the energies of individual spinon and holon can be 
obtained from GSEs of odd-L periodic chains with no holes 
and one hole, respectively. The even-L periodic chain with 
one hole contains the spinon-holon pair. A consistent finite- 
size definition of A^ for even L is: 
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where and 1 refer to the number of holes. It can be used 
directly to obtain a sequence of binding energies Al for dif- 
ferent L, which can then be extrapolated to L = oo limit. We 
will call this approach "method B". An alternative and sub- 
stantially more precise way to estimate Aoo is: (i) to obtain 
E X (L) by subtracting the extensive part of energy e a L from 
corresponding GSEs (e a is the energy per site of an infinite 
XXZ chain, known from Bethe-Ansatz (BA) [8[]), (ii) extrap- 
olate E X (L) to the infinite system size, and (iii) use Eq. ([3]). 
We will refer to this approach as "method A". In this method 
we use the exact L = oo value for the spinon energy E s that 
is also known from BA thus completely eliminating one 
of the sources of the finite-size effects. 

Method A. One of the problems is the "staggered" behavior 
of the holon and pair energies vs L, meaning that the data must 
be separated into two branches, defined by whether L for the 
pair (L — 1 for the holon) is divisible by 4. We refer to them 
as "4-even" and "4-odd" branches, respectively (see also Ref. 



2). Within each branch energy depends on the size smoothly. 
We found for both holons and pairs that although one of the 
branches may exhibit a non-uniform dependence on system 
size, the other one is always well-behaved, see Fig. [4] Thus, 
even though we have fewer points in the "good" branch, the 
quality of extrapolation using it is vastly improved. 

In fact, for the holon data, we can further advance this suc- 
cess by doubling the number of points in each branch using 
the "twisted" boundary conditions (BCs) |7[]. The GS of a 
holon in an infinite chain is degenerate for k = and k = ir. 
For a finite L, this degeneracy is lifted: in 4-even chains Eq is 
shifted up and E n is higher in the 4-odd ones (see Fig. [2]). The 
staggering occurs because the fc-space of odd-L chains does 
not contain the k = it point. This can be mitigated by flipping 
the sign of the hopping integral t, which shifts the momentum 
by it. Thus, we can construct the GSEs for both branches in 
every odd-L chain by combining the data as follows: 
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4-even branch. (5) 



The complementary data goes into 4-odd branch. An exam- 
ple of such reconstruction, along with the extrapolation to the 
infinite size, is shown in Fig. Eta). 

For the spinon-holon pair data, the separation into two 
branches also allows for a smooth E vs 1/L scaling. Further 
improvement using twisted BCs can be achieved by multiply- 
ing t with an a priori unknown L-dependent complex phase 
factor [7]. While such BCs can be easily handled by ED, 
our DMRG code requires extensive modifications to support 
them. Results featuring phase-adjusted data will be presented 
elsewhere fioll . 

Having significantly improved the scaling quality of the 
available GSE data, we need to choose the functional form 
for the E vs 1/L fit. We find that at not too large a GSEs ap- 
proach Eoo exponentially fast. Notably, similar exponential 
convergence of the energy of the finite XXZ chain at large 
L was derived by Woynarovich and de Vega (WdV) using BA 
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FIG. 4: The 4-even (circles) and 4-odd (squares) branches for (a) 
holon and (b) pair energies vs 1/L. J z — 4.0 and a — 0.4 (a = 
0.25, insets). Dashed curves are the fits described in text. 
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FIG. 5: Al data vs 1/L in method B, J z = 4.0 and a = 0.4. 
Branches are defined in Table U Dashed lines are guides to the eye, 
solid line is the 4-th order polynomial fit of the B2 branch. 



1 1 111 . Thus, we attempt to fit both holon and pair data with 
El = + e~ XL P n (l/L), where P n {x) is a polynomial of 
order n. Remarkably, the values of A we found for holons are 
in close agreement with WdV ones, so we used them in our 
holon fits. For pairs there is no such correlation, so we retain 
A's as free parameters. Example of pair energy extrapolation 
is shown in Fig. 15b). 

Finally, we are able to obtain high-precision values for the 
L = oo binding energies in the range < a < 0.5 using OJ. 
The data for J z = 4.0 are shown in Fig. [3] By studying the 
quality of the fits and the variation of Aoo depending on the 
fit type, we estimate the error to be less than 3% for a = 0.5 
and negligible for smaller values of a. For a > 0.5 parameter 
A drops below 1/L max , so the corresponding fits become less 
reliable. For a — 0.6 the error is of the order of 10%, and for 
a = 0.7 it exceeds 100%. In principle, our results at larger 
a can be improved by increasing L while keeping the DMRG 
precision the same (~ 10~H). However, as one can see from 
the inset in Fig. [3] the binding energy itself diminishes rapidly 
and it is likely to hit the DMRG accuracy around a ~ 0.8. 

Method B. The "direct" extrapolation of using Eq. 
(method B), does not rely on prior knowledge of BA results. 
Thus, it provides an important validity test for the results of 
method A. The A^ data is, again, staggered, and has to be 
separated into the 4-even and 4-odd branches. We also have a 
freedom in choosing the sign of t for the holon energies E\_ x 

L+l- 



and El , n . As a result, data splits into four different branches, 
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TABLE I: Subdivision of the numerical data for A(L) into different 
branches due to finite size effects in method B. 



defined in Table U An example of size dependence for dif- 
ferent branches is shown in Fig. [5] The best choice, corre- 
sponding to the "good" holon branch in method A, is branch 
B2. The results of extrapolation of B2 data using polyno- 
mial fit of maximum order are shown in Fig. [3] They are 
in close agreement with method A. Inspection of the terms in 
reveals that most significant finite-size effects in method B 
come from the spinon energies E^ ±i (eliminated in method 
A). Even though method B serves as an important validity 
test, its relative error is always larger than that of method A. 

Theory. Generally, the binding energy in ID should scale as 
—V 2 m, where V is the interaction and m is a particle mass. 
For the Ising limit this estimate gives A ~ — J%/t in agree- 
ment with the exact answer ([T). To extend our approach to 
the finite a, we use BSE formalism 11 1 211 in which two parti- 
cles with dispersions efe and uj q , and interaction V g>q ' create 
a bound state if their scattering amplitude has a pole. At the 
pole, BSE can be simplified to the Schrodinger equation for 
the pair wavefunction x(l) : 



x(q) 



1 



E - 
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(6) 



where P is the total momentum of the pair, and the pair energy 
E = A + e + oj is the binding energy A relative to the band 
minima of the the particles eo = mm[efe] and w = min[ui q }. 
In the Ising limit, = — 2t cos k, u> q = J z /2, and V q . q ' = —ojq, 
so Eq. © yields a dispersionless (P-independent) bound state 
with A given by (HJ. 

At a > spinon becomes mobile and the interaction V 
and holon dispersion change. The holon mass renormaliza- 
tion was found to be insignificant throughout the anisotropic 
regime < a < 1 [7]. On the other hand, spinon dis- 
persion changes drastically as it evolves from the gapped, 
immobile excitation with energy uj q = J z /2 at a = to 
relativistic, gapless excitation with ui q = J z (ir/2)\ cosq\ in 
the isotropic limit. The spinon dispersion for XX Z model, 
shown by solid lines in Fig. |2ja), is known exactly from 
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BA: Lo a = c\J 1 — k 2 sin 2 q, for parameters c and k see 
Ref. |9l While c/J z changes almost linearly between 1/2 
and 7r/2 as a goes from to 1, k varies steeply from 
to 1, such that the spinon gap uj s = u} m in becomes small 
already for a > 0.5 and approaches zero exponentially: 

lo s k, 4cexp ^— tt 2 yJa/%{\ — a)J as a — > 1. This makes 
the spinon spectrum "quasi-relativistic" in this regime, w g = 
J icq) 2 + uj 2 . The spinon also becomes very light, its mass 
going from m s ^a -1 at a <C 1 to to s ~cj s at a >0.5, see Fig. 
|2ja). Such a transformation of the spinon spectrum strongly 
affects the spinon-holon binding. Note that, when the spinon 
is much lighter than the holon, m s <C m/, ~ (2t) _1 , the role 
of the latter in pairing must be secondary. 

The remaining question is the spinon-holon interaction. 
The binding problem can be solved rigorously to order 0(a). 
In the small-a regime, changes to the holon and the AF 
GSEs are of order 0(a 2 ), while the spinon energy changes 
in the order 0(a): uj q = uj° + 8u q , where u° = J z /2 
and 5u! q = aJ z cos2q. At a > spinon GS momentum 
is ±7r/2, so the bound state should also have a finite mo- 
mentum P — ±7r/2, in agreement with the numerical data. 
Since the energy is lowered when the AF domain walls as- 
sociated with spinon and holon cancel each other, the inter- 
action can be written as a "contact" attraction of strength 
V° = —J z /2. This leads to a relation between interaction 
in the momentum space and spinon dispersion, valid to order 
0(a): Vq t q> = — u — (Sujq + SLJqi)/2. Using this interaction, 
spinon energy as above, and the "bare" holon energy e? in 
Eq. (O yields: A = Ao(l — Act), with Ao given by (Q~|i and 
A\ A | = + J z - The slope A varies from 4 to 2 for 

< J z < 00. This 0(a) result for A at J z = 4.0 is shown in 
Fig. [5]by the dashed line. It is in excellent agreement with the 
numerical data in the small-a regime. 

Based on the above analyzis, we propose the general form 
of the spinon-holon interaction in momentum space: V Qiq r = 
~^jui q LOqi . Using this Vq >q ', spinon energy from BA, and 
"bare" holon dispersion, Eq. (O is transformed into: 



g A - (e q+P - e ) - (uj q - u) ) ' 

Solving this equation gives the dependence of the binding en- 
ergy A on anisotropy a shown as a solid line in Fig. [3] Not 
only does this equation yield our small-a results, but it also 
provides a very close agreement with the numerical data for 
all the values of J z and for all a we can access numerically. 
This agreement makes the validity of our spinon-holon inter- 
action ansatz very plausible. 

As the binding energy becomes small at a — > 1, pair- 
ing is determined by the long-wavelength features of the 
dispersions and interaction. Within the qualitative picture 
of pairing in ID, both the characteristic low-energy interac- 
tion and the spinon mass become proportional to the spinon 



gap U! s , Thus, one expects A ~ —V 2 m ~ — ui^ ~ 
— exp ^— 37r 2 / y/8(l — on . We can derive this behavior 
from Eq. (0 explicitly: A w —T>(J Z , a)wf , where the expo- 
nential behavior is determined solely by the spinon gap. The 
holon energy scale is secondary as it only enters the "regu- 
lar" pre-factor T>(J Z , a). Altogether, this explains the quick 
(exponential) drop-off of A at intermediate a. 

One can see from the asymptotic expression that the bind- 
ing energy vanishes in the isotropic limit together with the 
spinon gap. Thus, our spinon-holon interaction ansatz pro- 
vides a natural and simple explanation of the non-zero binding 
but no bound state at a = 1: it is possible because the inter- 
action of the holon with the long-wavelength spinon vanishes 
together with the spinon energy. Then the pairing is not strong 
enough to produce a bound state. 

Conclusions. To summarize, we have studied the spinon- 
holon interaction in the anisotropic t-J model. We have 
demonstrated that the finite-size extrapolation of the ED and 
DMRG data of very high accuracy is possible and it can pro- 
vide reliable values of the spin-holon binding energy. These 
data are in excellent agreement with the theory based on the 
BSE with the contact spinon-holon interaction. The theory 
provides a coherent, simple, and consistent explanation of the 
binding for all values of anisotropy. Within the offered pic- 
ture, the holon is attracted to the "slow" spinon in the a — > 
limit while at a > 0.5 the "fast" spinon can only be bound 
weakly. We also offer an explanation of the non-zero attrac- 
tion but no bound state in the isotropic t-J model. 
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